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(Dated: February 7, 2008) 

We report on a search for the pair production of scalar leptoquarks, LQ, using 191 pb~^ of proton- 
antiproton collision data recorded by the CDF experiment during Run II of the Tevatron. The 
leptoquarks are sought via their decay into a neutrino and quark yielding missing transverse energy 
and several jets of large transverse energy. No evidence for leptoquark production is observed, and 
limits are set on cr{pp LQLQX vvqqX). Using a next-to- leading order theoretical prediction 
of the cross section for scalar leptoquark production, we exclude first-generation leptoquarks in the 
mass interval 78 to 117 GeV/c^ at the 95% confidence level for BR(L(5 ^ vq) = 100%. 

PACS numbers: 12.60.-1, 13.85.Rm, 14.80.-j 



The remarkable symmetry between quarks and leptons 
in the standard model (SM) suggests that some more 
fundamental theory may exist, which allows interactions 
between them. Such interactions are mediated by a new 
type of particle, a leptoquark |1|], which carries both lep- 



ton and baryon number. A leptoquark is a color-triplet 
boson with spin or 1, and has fractional electric charge. 
Leptoquarks are predicted in many extensions of the SM 
(e.g. grand unification, technicolor, and supersymmetry 
with i?-parity violation). The Yukawa coupling of the 
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leptoquark to a lepton and quark and the branching ratio 
to a charged lepton, denoted by /?, are model dependent. 
Usually it is assumed that leptoquarks couple to only one 
generation to accommodate experimental constraints on 
flavor-changing neutral currents Q , which allows one to 
classify leptoquarks as first-, second-, or third-generation. 
In pp collisions, leptoquarks can be produced in pairs 
via the strong interaction through gg fusion or qq anni- 
hilation. The production rate for scalar leptoquarks is 
essentially model-independent and is determined by the 
known QCD couplings and leptoquark mass. 

We report on a search for pair production of scalar 
leptoquarks, with LQ decaying to vq, resulting in a jets 
and missing transverse energy (I^t) topology. We use 
191 ±11 pb~^ 3 of pp collision data at a center-of-mass 
energy of 1.96 TeV recorded by the Collider Detector at 
Fermilab (CDF) during the Tevatron Run II. This analy- 
sis is sensitive to leptoquarks of all three generations with 
/3 « 0. The previous lower mass limit of 98 GeV/c^ 
on first-generation leptoquarks in this final state was set 
by the D0 Collaboration. The CDF Collaboration has 
also published lower mass limits of 123 GeV/c^ and 
148 GeV/c^ respectively on second- and third-generation 
leptoquarks in the I^t plus heavy-flavor jets final state. 
Limits on leptoquark production from the Tevatron Run 
I and HERA experiments as of 1999 are summarized in 
i 

CDF is a generalj)urpose detector that is described 
in detail elsewhere [g]- The components relevant to this 
analysis are briefly described here. The charged-particle 
tracking system is closest to the beam pipe, and consists 
of multi-layer silicon detectors and a large open-cell drift 
chamber covering the pseudorapidity Q region < 1. 
The tracking system is enclosed in a superconducting 
solenoid, which in turn is surrounded by a calorimeter. 
The CDF calorimeter system is organized into electro- 
magnetic and hadronic sections segmented in projective 
tower geometry, and covers the region \r]\ < 3.6. The elec- 
tromagnetic calorimeters utilize a lead-scintillator sam- 
pling technique, whereas the hadron calorimeters use 
iron-scintillator technology. The central muon-detection 
system, used for this analysis, is located outside of the 
calorimeter and covers the range \ri\ < 1. 

This search centers on selecting events with large I^t 
and a pair of jets that are acoUinear in the transverse 
plane, because of the neutrinos in the final state. The 
I^T y is defined as the energy imbalance in the plane 
transverse to the beam direction. A jet is defined as 
a localized energy deposition in the calorimeter and is 
reconstructed using a cone algorithm with fixed radius 
AR = ^/ArfT'Acj^ ^ OA in rj - <j) space We cor- 
rect jet Et measurements and I^t for detector effects. 

The data sample for this analysis was collected using 
an inclusive I^t trigger, which is distributed across three 
levels of online event selection. In the first and second 
levels of the trigger, I^t is required to be greater than 



25 GeV and is calculated by summing over calorimeter 
trigger towers (IQ] with transverse energies above 1 GeV. 
At Level-3 I^t is required to be greater than 45 GeV and 
is recalculated using full calorimeter segmentation with a 
tower energy threshold of 100 MeV. We use events from 
the inclusive high-p^ lepton (e or /i) samples to measure 
the trigger efficiency directly from data. To reduce sys- 
tematic effects associated with the online trigger thresh- 
old, we select events offline with I^t > 60 GeV, where 
the trigger is fully efficient. 

The event electromagnetic fraction (Fem) and charged 
fraction (Fch) 01 nsed to remove events associated 
with beam halo and cosmic ray sources. We reject events 
that contain little energy in the electromagnetic section 
of the calorimeter or that have mostly neutral-particle 
jets, by requiring F^m > 0.1 and Fch > 0.1. There are 
148,462 events in our analysis sample after the initial 
selection. 

The dominant backgrounds to the leptoquark search in 
the jets and I^t signature are QCD multi-jet production, 
W and Z boson production in association with one or 
more jets, and top quark pair production. The ALPGEN 
generator 0| was used for the simulation of the W and 
Z boson plus parton production, with HERWIG used 
to model parton showers. We use the exclusive Z{^ 
ee)-|-l jet sample to determine a scale factor between data 
and simulation, and apply this factor to all W/Z+jets 
simulation samples, herwig was also used to estimate 
the contribution from tt production. 

Data selection requirements were chosen to maximize 
the statistical significance of the leptoquark signal over 
background events based on studies of simulated event 
samples before the signal region data were examined. In 
addition to I^t > 60 GeV, the signal region is defined by 
requiring that the two highest Et jets (i?;^^ > 40 GeV, 
Ei^i^ > 25 GeV) be in the central region |r/| < 1. A third 
jet with Et > 15 GeV and |ry| < 2.5 is allowed, and we 
veto events with any additional jets with Et > 15 GeV 
and 1 77 1 < 3.6. To reject events with I^t resulting from 
jet energy mismeasurement, we require that the opening 
angle in the transverse plane between the two highest Et 
jets satisfy 80° < A(?!)(ji, ja) < 165°. The ^t direction 
must not be parallel to any of the jets; we require the 
minimum azimuthal separation between the direction of 
the jets and ^t to satisfy 30° < min A(/)(j, ^r) < 135°. 
The I^T also must not be antiparallel to the leading Et 
jet: 100° < A(/)(ji,^T) < 165°. These criteria reject 
most of the QCD multi-jet background events. To re- 
duce the background contribution from W/Z +}ets and 
tt production, we reject events with one or more iden- 
tified leptons with Et > 10 GeV (electron candidates) 
or pt > 10 GeV/c (muon candidates). Criteria similar 
to those in 0| are used to identify the leptons. To fur- 
ther reduce this background we require each jet not to 
be highly electromagnetic (jet electromagnetic fraction 
< 0.9) and to have 4 or more associated tracks for cen- 
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tral jets (I77I < 1). 

Two methods are employed to estimate the QCD 
multi-jet contribution in the signal region directly from 
the inclusive I^t data sample. Among all the offline anal- 
ysis selection requirements, the azimuthal angular sepa- 
ration requirement between the I^t direction and a jet is 
most effective at removing QCD multi-jet events. There- 
fore, for the first method, in addition to the signal re- 
gion we define a region which is rich in QCD multi-jet 
events by requiring that a jet is close to the I/^t direc- 
tion (20° < minA(^(j,^T) < 27°). Studies of simulated 
QCD multi-jet samples show that the shape of the I^t 
distribution in this region is similar to the I^t distribu- 
tion in the signal region. We use Ifx and min A(/)(j. I/Jt) 
requirements to define four kinematic regions: 

A) 50 < < 57 GeV, 20° < min A0(j, ^t) < 27°. 

B) > 60 GeV, 20° < min A(/)(j, ^t) < 27°. 

C) 50 < < 57 GeV, 30° < min A</>(j, ^t) < 135°. 

D) > 60 GeV, 30° < minA0(j,^T) < 135°. 

The regions A, B and C are used to extrapolate the 
QCD multi-jet contribution into the signal region D: 
Nd = ^Nq, where N^, Nb, and Nc are the remain- 
ing number of events in regions A, B, and C, after the 
W/Z+iets and tt contributions have been subtracted. 
For the second method, the combined selection require- 
ment efficiency is measured as a function of I^t in an 
independent inclusive jet sample at low I^t ■ The ex- 
trapolated results of this measurement is then applied to 
the inclusive I^t sample after the W/Z+jets and ti con- 
tributions have been subtracted. We predict 15.0 ± 8.0 
and 21.5 ± 12.4 multi-jet events for the first and second 
methods respectively. We take the weighted average and 
uncertainty of the two methods as our estimate of the 
multi-jet background. 

We check the simulation predictions for W/Z+jets 
with data in a control region, which is defined by re- 
quiring, in addition to 2 or 3 jets, I^t > 60 GeV and 
at least one electron or muon. We observe 144 events in 
our inclusive I^t sample, which is in excellent agreement 
with 154.3 ± 27.9 events predicted from SM processes. 

The total detection efficiency (clqi) for the first- 
generation scalar leptoquark (LQi) signal is estimated 
using the PYTHIA event generator [ij , and the CDF de- 
tector simulation program. The PYTHiA underl ying event 
simulation was tuned to reproduce CDF data |lfil|. The 
samples were generated using the CTEQ5L [13I parton 
distribution functions (PDF), with the renormalization 
and factorization scales set to /i = m^Qj. Table lists 
the total detection efficiency elQi and the corresponding 
total fractional uncertainty Stot for various lep toq uark 
masses. Also listed are the NLO cross sections [ig cal- 
culated for two choices of the fj, scale. The systematic 



TABLE L Summary of the first-generation scalar leptoquark 
detection efficiency (elQi), the relative uncertainty on de- 
tection efficiency (Jtot), and the next-to-leading order cross 
section (ctnlg) for two choices of the renormalization scale as 
functions of leptoquark mass. 



o-NLo(pb) 

rriLQi (GeV/c^) e^Qi ^tot (%) fJ, ^ tulQi ^ = 2mz,Qi 



75 


0.0073 


29 


69.4 


58.8 


80 


0.0113 


26 


49.2 


41.5 


90 


0.0187 


23 


26.0 


22.1 


100 


0.0300 


20 


14.6 


12.5 


110 


0.0431 


16 


8.4 


7.4 


115 


0.0482 


15 


6.7 


5.8 


125 


0.0590 


15 


4.2 


3.6 


150 


0.0828 


13 


1.4 


1.3 


175 


0.1010 


12 


0.57 


0.51 



uncertainty on the signal acceptance includes the uncer- 
tainties due to modeling gluon radiation from the initial- 
state or final-state partons (10%), and the choice of the 
PDF (4%). The limited size of the leptoquark simula- 
tion samples gives a 3% statistical uncertainty. The sig- 
nal acceptance uncertainty due to the jet energy scale 
varies from 4% to 26%, and the uncertainty on the lu- 
minosity is 6%. The uncertainty on the trigger efficiency 
is 1%. The theoretical uncertainties on the renormaliza- 
tion and factorization scales are not included here, since 
we conservatively choose the NLO cross section setting 
fi — 2mLQi to extract the limits on leptoquark mass. 
This choice of scale is found to reduce the cross section 
prediction by 15% relative to /i = rriLQ^ [l^ . 

In the signal region, we expect 118.5 ± 14.5 events 
from SM processes and observe 124 events. The pre- 
dicted backgrounds from SM processes are summarized 
in Table ^1 In Figure ^ the predicted Ifx distribution 

TABLE II: The number of expected events from various SM 
sources in the leptoquark signal region. The first uncertainty 
is from the limited simulation statistics and the second is from 
the various systematics. 



Source Events expected 

W{^eu)+]ets. 6.1 ±1.4 ±1.5 

W{^^.v)+]eis 21.7 ±2.3 ±2.8 

W{^Tv)+]ets 28.4 ±3.8 ±4.1 

Hn)+]ets 1.1 ±0.2 ±0.2 

Z(^rr)+jets 0.9 ± 0.2 ± 0.2 

zl-^vv)+]els 39.1 ±2.8 ±3.6 

tt 4.3 ±0.4 ±0.3 

QCD 16.9 ±6.7 

Total Events 118.5 ± 14.5 



is compared with the distribution observed in data. No 
evidence for leptoquark production is observed. We cal- 
culate the upper limit at the 95% confidence level (C.L.) 
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FIG. 1: The Ifr distribution in the leptoquark signal re- 
gion for data (solid points) compared to SM background 
(shaded histograms). Also shown is the expected distribu- 
tion arising from leptoquark production and decay at a mass 
of 115 GeV/c^ (hatched histogram). 



on the pair production cross section times the square of 
the branching ratio of the leptoquark to a quark and a 
neutrino using a Bayesian approach pjj with a fiat prior 
for the signal cross section and Gaussian priors for ac- 
ceptance and background uncertainties. The upper limit 
on the cross section times (1-/3)^ is shown in Figure El 
and is compared with the theoretical cross sections. The 
theoretical cross sections for scalar leptoquark produc- 
tion have been calculated at NLO using CTEQ5M flTj 
PDFs. 

In conclusion, we performed a search for leptoquarks 
in the jets and Ifr topology using 191 pb~^ of CDF Run 
II data. No evidence for leptoquarks is observed. We 
set an upper limit on the production cross section at the 
95% C.L. Assuming a leptoquark decays into a neutrino 
and quark with 100% branching ratio, we exclude the 
mass interval from 78 to 117 GeV/c^ for first-generation 
scalar leptoquarks. This extends the previous limit for 
the first-generation scalar leptoquark of 98 GeV/c^ 0. 
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FIG. 2: The upper limit on the cross section times squared 
branching ratio for scalar leptoquark production in the jets 
and Ifr topology. Also shown is the NLO cross section for 
/3 = for 2 choices of the factorization/renormalization scale: 
M = mLQi, /i = 2 mLQi. 
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